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1.  Introduction: 


During  the  contract  period,  the  principal  investigator  and  his  graduate  associates  Mr. 
C.Y.Ching  and  Ms.  B.  Ayotte  worked  on  several  problems  related  to  lead-induced  convection,  and 
oceanic  convection  in  general.  The  results  of  these  studies  are  outlined  below.  Mr.  Ching  is  in  the 
process  of  Hnishing  his  Ph.D.,  and  his  thesis  will  be  available  as  a  report.  Other  relevant 
publications  originated  from  the  P.L's  group  during  the  contract  period  are  listed  at  the  end  of  this 
report 


2.  Research  on  Oceanic  Convection 

Because  of  large  space  and  time  scales  involved,  oceanic  convective  motions  (which  may  take 
the  form  of  plumes  and  thermals)  are  expected  to  be  influenced  by  the  background  rotation.  These 
rotational  effects  may  arise  either  due  to  the  earth's  rotation  or  due  to  the  presence  of  a  parent 
cyclone.  Thus,  the  understanding  of  oceanic  convection  requires  a  good  knowledge  on  the 
influence  of  rotation  on  convection  from  isolated  sources. 

(a)  Plumes 

If  a  point  plume  is  released  from  a  source  placed  within  a  fluid,  in  the  absence  of  rotation,  it  wiU 
first  accelerate  to  a  maximum  speed  and  then  decelerate.  If  the  source  diameter  of  the  plume  is  do, 

then  it  is  possible  to  argue  that  its  evolution  will  be  governed  by  the  total  buoyancy  defined  as  B  = 
ndo^PoUo/^’  where  Po  is  the  initial  buoyancy,  uq  is  the  initial  velocity  and  do  is  the  source 


diameter  of  the  plume.  Thus,  it  is  possible  to  write  the  r.m.s.  velocity  u  and  the  integral 
lengthscale  t  of  turbulence  within  the  plume,  at  a  distance  z  away  from  the  source,  as 


(1) 

(2) 


and 

where  fi,  f2  . . .  are  functions  and  aj,  a2  ...  are  constants.  For  point  plumes  do/z  0,  and  it 
is  possible  to  write 
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Further,  the  depth  h  to  which  the  plume  descends  at  a  time  t  after  the  release  can  be  written  as 

h  =  (Bt3)  l%3(do/(Bt3)lM),  (5) 

where  purely  dimensional  arguments  have  been  used.  For  point  plumes  (5)  becomes 

h  =  a3(Bt3)l/4  (6) 

When  the  plume  descends  in  a  rotating  fluid,  it  is  expected  that  it  will  be  influenced  by  Coriolis 
and  centrifugal  forces  some  time  after  the  onset  of  the  motion.  Since  the  mean  velocity  of  the 
plume  is  vertical  (in  the  direction  of  the  rotational  vector),  the  rotational  influence  will  act  mainly  on 
the  plume  turbulence  rather  than  the  mean  vertical  flow.  Previous  works  of  the  P.I.  have  shown 
that  the  turbulence  is  affected  by  the  background  rotation  when  the  inertia  forces  are  of  the  same 
order  as  the  Coriolis  forces,  viz., 

»^»2Qu,  (7) 

or,  when  the  turbulent  Rossby  number  drops  below  a  critical  value, 

Ro  =  4L  <  Rocr  =  04  ,  (8) 

Cht 

where  Cl  is  the  rate  of  rotation.  Using  (3),  (4),  and  (8)  for  point  plumes,  the  critical  depth 
where  the  rotational  effects  inhibit  the  growth  of  the  pluizx;  can  be  wrinen  as 

Since  the  width  of  the  plume  b  at  z  <  can  be  vmtten  as 

b  =  zf4(%),  (10) 


or  for  point  plumes 


b  =  Osz, 


the  critical  width  of  a  point  plume  be  at  the  onset  of  rotational  effects  is 
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Using  (6)  and  (9),  the  time  tc  for  the  onset  of  the  rotational  effects  can  be  estimated  as 

fltc=a7. 
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In  oceanic  cases,  the  descending  plumes  eventually  intersect  with  a  horizontal  boundary  (say 
located  at  a  distance  H  below  the  source),  and  it  is  instructive  to  consider  the  evolution  of  the 
plume  once  it  impinges  on  the  boundary  and  spreads  horizontally  as  a  gravity  current.  Two  cases 
should  be  considered. 

(a)  The  plume  is  not  influenced  by  the  rotation  prio*  to  the  arrival  at  the  boundary  (zc  »  H);  In 
this  case,  the  plume  parameters  at  the  impingement  can  be  calculated  using  the  governing  variables 
B,  H,  and  do-  For  the  simple  case  of  a  point  plume,  the  buoyancy  of  the  plume  at  the  interface  |3 
can  be  written  as 


or  for  a  plume  with  a  fmite  diameter 


(14) 


(15) 


Upon  impingement,  the  resulting  gravity  current  flows  horizontally  until  the  baroclinic  instabilities 
are  developed.  The  velocity  of  the  gravity-current  nose  can  be  evaluated  by  a  balance  between 
horizontal  pressure  gradient  forces  and  inertia  forces  as 

c-fehg)*'^  ,  (16) 

where  hg  is  the  height  of  the  gravity  current.  On  dimensional  grounds,  for  a  point  plume, 

hg~H  (17) 


and  hence 

c~(Ph)^'^  .  (18) 

If  one  assumes  that  the  gravity  current  becomes  unstable  and  breaks  into  a  number  of  eddies, 
whence  it  propagates  to  a  radius  rc  proportional  to  the  Rossby  radius  of  deformation  c/O,  it  is 
possible  to  write 

where  ag  is  a  constant.  On  the  other  hand,  if  a  plume  with  a  finite  source  diameter  do  is 
considered,  the  dimensional  arguments  yield 
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(20) 
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However,  in  the  limit  do  0,  this  should  yield  (19)  and  hence  it  is  possible  to  argue  that  f6 
should  be  independent  of  H/(B/Q3)1/4^  viz., 

_  fH.I 

(21) 
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£b)  zc  <  H:  Since  the  plume  growth  is  arrested  by  the  rotation  at  zc,  it  is  possible  to  assume  that 
the  buoyancy  at  z  =:  H  is  the  same  as  that  at  z  =  zc-  Using  an  expression  equivalent  to  (14),  it  is 
possible  to  write 

(22) 

and  equivalent  to  (18), 

c~(pZcr-(BQf'\  (23) 

Thus  the  radius  at  which  the  flow  (gravity  current)  becomes  unstable  is  given  by 


rc» 


(24) 


This  suggests  the  possibility  of  the  gravity  current  becoming  unstable  from  the  onset, 
(b)  Thermals 


For  thermals  released  from  a  point  source,  it  can  be  shown  that  the  governing  parameter  would 
be  the  total  initial  buoyancy  qj.  Previous  self-similarity  arguments  and  laboratory  experiments 

suggest  that,  in  the  absence  of  background  rotation,  the  vertical  distance  traveled  by  the  thermal  h 
and  its  width  b  at  any  time  t  during  the  decelerating  phase,  should  satisfy 

h  =  ciqx^/^^^.  (25) 

b  »c2h ,  (26) 

where  ci,  C2  are  constants.  Similarly,  the  turbulent  quantities  pertinent  to  a  descending  thermal 
can  be  written  as 


Oji/z 

z  ’ 

(27) 

and 

6-z. 

(28) 
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It  is  expected  that  the  thennals  feel  the  rotational  effects  when  (8)  is  satisfied  or  when 


and 


b^«C4p^ 


(29) 

(30) 


Upon  the  thermal  reaching  the  horizontal  surface,  it  is  expected  to  spread  as  a 
become  unstable  to  shed  baroclinic  eddies.  As  for  the  case  of  plumes, 
considered. 

(a)  zc  »  H:  For  the  simple  case  of  a  point  thermal,  the  buoyancy  of 
horizontal  boundary  can  be  estimated  as 


gravity  current  and 
two  cases  can  be 

the  thermal  at  the 


(31) 


and  the  frontal  propagation  velocity  is  given  by 


(32) 


Since  the  characteristic  height  of  the  thermal  at  any  frontal  radius  r  can  be  estimated  as  r^hg^H^,  it 
is  possible  to  write 

or  r2  -  r?  »  qpt,  (34) 

where  r©  is  the  impact  radius  of  the  thermal  and  t  is  the  time  elapsed  after  the  impingement.  The 


radius  at  which  the  gravity  current  becomes  unstable  can  be  evaluated  as  c/Q  or 

>a' 


1 


(35) 


(b)  Zc  «  H:  In  this  case  the  thermal  growth  is  arrested  by  the  rotation  at  zc  and  it  can  be 
assumed  that  the  buoyancy  at  z  =  H  is  the  same  as  that  at  z  =  zq.  The  buoyancy  of  the  thermal  at 
the  impingement  is  given  by  (qT^^)^/^,  and  thus 

(36) 


Laboratory  experiments  were  performed  by  releasing  either  a  turbulent  plume  or  a  thermal 
(made  out  of  salty  water)  to  a  homogeneous  (or  a  two-layer)  fluid  contained  in  a  rotating  tank. 
Standard  flow  visualization  and  particle-tracking  systems  were  used  to  obtain  information  on  the 
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flow  evolution.  Plumes  of  different  diameters  were  also  considered  to  delineate  the  effects  of  the 
source  diameter.  In  the  former  case,  the  starting  plume  sank  into  the  underlying,  less-dense, 
homogeneous  layer  of  water  forming  a  growing  3-D  mrbulent  layer.  When  the  turbulent  front 
reached  a  transition  depth  zq,  the  rotational  effects  became  dominant  and  the  lateral  growth  of  the 
plume  was  found  to  be  inhibited.  Ruid  motions  appeared  to  be  2-dimensional  in  planes 
perpendicular  to  the  axis  of  rotation.  At  this  time,  the  descent  of  the  plume  was  somewhat 
decelerated  and  2-D  vortex  structures  were  found  to  form  beneath  the  3-D  turbulent  mixed  layer. 
These  vortices  then  penetrated  downward  to  produce  vtHiex  columns.  Once  the  plume  reaches  the 
bottom  of  the  tank,  it  spread  slowly  along  the  bottom  forming  2-D  vortex  columns.  Figure  1 
shows  a  plot  of  the  frontal  position  versus  time,  using  the  non-dimensionalized  parameters 
h/(B/0^)l/^  and  Ot.  Note  the  initial  descent  with  3/4  slope.  Then  when  h/(B/03)^/^  is 
approximately  3.3  at  Ot=  2.3,  the  slope  decreases  showing  a  change  in  speed  of  the  starting  plume 
cap.  The  corresponding  variation  of  the  width  is  shown  in  Figure  2.  Note  the  complete  arrest  of 
the  horizontal  growth  of  the  plume  at  a  time  Ch  »S.S.  The  calculated  parameters  based  on  the 
experiments  are  aze  *  3.3  and  abc  *  1*8. 

When  the  result  zc  =  3.3(B/Q3)1/4  Js  applied  to  the  ocean  with  typical  values  of  B  = 
412  m^s*3  and  Q  =  7.14  x  lO'^  rads'^,  it  is  possible  to  get  zc  -19  km.  This  is  obviously  greater 
than  the  depth  of  the  ocean,  and  hence  the  geophysically  interesting  case  would  be  H  <  zc  .  The 
observations  made  during  the  experiments  with  H  <  z^  show  that  the  plume  impinges  on  the 
bottom  of  the  tank  and  propagates  horizontally  as  a  gravity  current,  until  it  grows  to  a  certain 
lengthscale,  re-  Thence  the  flow  becomes  baroclinically  unstable  and  breaks  up  into  isolated 
eddies.  The  critical  radius  rc  at  which  the  eddies  are  shed  was  measured  and  was  plotted  as  a 
function  of  (B/fHfl^))^^;  a  forced  fit  line  of  slope  one  indicated  ®g  =  1.3. 

Some  experiments  were  carried  out  in  non-rotating  fluids  to  investigate  the  effect  of  the 
source  diameter  on  plume  properties.  According  to  (5),  h  =  (Bt3)l/^3(d<^(Bt3)l/^ );  and  when  do 
->0  or  do/(Bt3)l/4  -»0  (large  t),  h  =  a3(Bt3)l/^.  On  the  other  hand,  when  the  source  is 
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homogeneous,  the  rate  of  descent  becomes  h  =  ai|Bt^/cl^)l/2  (second  kind  of  self  similanty  when 
do  Thus  we  expect  the  characteristic  descent  law  of  the  starting  plume  to  change  from 

to  t3/2  as  the  source  diameter  changes  from  do~^  0  to  do”^®**-  In  other  words,  a  plot  of  h/(Bt3)l/4 
should  change  from  a  constant  aj  in  the  limit  do~^  0  to  a*(Bt^/d^)  1/4,  when  do  is  increased. 

Experiments  were  carried  out  with  varying  initial  source  diameters  do,  and  the  results  were  used 
to  evaluate  the  condition  under  which  the  initial  plume  diameter  can  be  neglected.  The  results  for 
different  do  values  initially  (at  small  t)  followed  t3/4,  but  then  clear  departures  could  be  seen. 
Accordingly,  the  point-plume  source  approximation  is  valid  when  the  plume  front  satisfies  h  > 
11.7do. 

Experiments  were  also  carried  out  to  investigate  the  effect  of  do  on  the  propagation  of 
gravity  currents  induced  by  the  impingement  of  plume  on  the  interface,  in  the  presence  of 
background  rotation;  as  was  discussed  earlier,  the  radius  at  which  the  gravity  current  becomes 
unstable  should  be  given  by  rc/(B/Hfl3)l/3  =  f7[do/H],  when  H  <  zc.  The  results  show  that, 

when  the  initial  diameter  is  unimportant,  the  non-dimensional  critical  radius  is  approximately 
1.4(B/Hfi3)l/3  When  do  becomes  importam  at  do/H  *  0.12,  the  slope  of  the  graph  showed  a 
decrease  according  to  rc/(B/(HQ^))^/3  =  -29(do/H)  +  5. 

3.  Research  on  lead-induced  Convection 

The  polar  ice  cap  often  cracks  and  forms  long,  narrow  channels  of  open  water  which  are 
known  as  leads.  These  openings  play  an  important  role  in  the  heat  budget  and  circulation  of  polar 
oceans.  In  the  winter  the  opening  of  a  lead  is  associated  with  the  onset  of  turbulent  convection 
above  and  below  the  lead,  because  of  the  refreezing  of  surface  water  and  rejection  of  heat  into  the 
atmosphere,  respectively.  Lead-induced  convection  is  associated  with  a  rich  variety  of  fluid- 
dynamical  phenomena,  yet  studying  them  in  field  situations  is  not  easy  due  to  the  unpredictability 
of  their  occurrence  and  short  resident  times.  As  such,  laboratory  experiments  and  numerical 
models  can  play  a  major  role  in  lead-related  studies.  The  purpose  of  this  study  was  to  review  the 
results  of  some  previous  laboratory  experiments,  which  appear  to  be  of  importance  in  guiding  and 
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interpreting  field  experiments,  and  to  carry  out  new  labcH^tory  experiments  dealing  with  lead- 
induced  motions. 

To  this  end,  laboratory  experiments  were  carried  out  to  investigate  the  interaction  between 
turbulent  line  buoyant  plumes  and  sharp  density  interfaces  with  the  hope  of  using  the  results  to 
interpret  oceanic  observations  pertinent  to  crack  openings  in  the  polar  ice  cap  (leads).  The 
experiments  were  carried  out  by  releasing  a  turbulent  line  plume  (with  buoyancy  flux  per  unit 
length  qo)  onto  the  surface  of  a  two-layer  fluid,  and  the  evolution  of  the  plume  was  nranitored. 
The  plumes  initially  descend  as  in  a  homogeneous  fluid,  impinge  on  the  density  interface,  and  then 
spread  horizontally  as  gravity  currents.  The  results  are  summarized  below. 

(i)  The  'cap'  of  the  starting  plume  initially  descends  as  in  a  non-stratified  fluid,  and  its 
position  h  and  width  b  at  a  given  time  t  is  given  by  h  =  Qo'^t  and  b  ~  0.66  h.  The  half  width  of 
the  plume  that  immediately  follows  the  cap  was  found  to  be  ^  =  0. 16  h. 

(ii)  The  evolution  following  the  impingement  of  the  plume  cap  on  the  density  interface 
was  found  to  be  governed  by  the  bulk  Richardson  number  Ri  =  Ab^iyqo2/3,  where  is  the  half 
width  i  of  the  plume  at  the  interface  and  Ab  is  the  buoyancy  jump  across  the  interface.  When  Ri  > 
10,  the  plume  cap  was  found  to  be  separated  into  two  hetrostrophic  vortices,  which  travel  a 
distance  of  x^.  =  (3  -  4)/d,  before  starting  their  gravitational  collapse.  The  resulting  gravity  current 
showed  some  similarities  to  those  propagating  over  solid  surfaces,  such  as  K-H  instabilities  behind 
the  head,  but  its  nose  angle  was  much  steeper.  When  Ri  >  5,  the  penetration  of  the  plume  into  the 
lower  layer  was  found  to  be  small,  of  the  order  of  the  thickness  of  the  interface.  When  1<  Ri  <  5, 
the  gravity  current  that  is  formed  by  the  plume  fluid  rebounding  off  the  distorted  interface  showed 
qualitative  differences  with  the  case  Ri  >  10.  The  vortex  nature  was  weak  within  the  current  and 
the  nose  angle  was  close  to  the  theoretical  value  predicted  by  the  inviscid  theory;  however,  no 
significant  mixing  was  present  at  the  bottom  boundary.  When  0.6  <  Ri  <  1,  mixing  was  prevalent 
between  the  gravity  current  and  the  surrounding  fluid  and  the  gravity  current  was  syrrunetric  with 
respect  to  the  interface.  At  much  lower  Ri  (<  0.5),  the  plume  was  found  to  penetrate  deep  into  the 
interface,  impact  on  the  tank  bottom,  and  spread  horizontally  while  rising  slowly. 
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(iii)  The  average  depth  5  of  penetration  of  the  piume  into  the  lower  heavy  layer  was  found 
to  be  given  by  =  4  Ri'l,  for  Ri  <  5.  At  large  Ri,  this  power  law  was  closer  to  Ri'*/2;  in  this 
case  5  is  smaller  than  the  interfacial  thickness,  and  hence  the  penetration  depth  follows  a  different 
power  law  behavior  than  that  of  a  two-layer  fluid. 

(iv)  The  outward-flowing  gravity  current  generated  by  the  plume  impingement  on  the 
interface  showed  a  constant  velocity  Uf  =  0.95  for  Ri  >  5.  At  low  Ri,  the  result  followed 
Uf/qo*^  =  0.7  Ri0  l7.  -Thg  average  thickness  of  the  current  is  given  approximately  by  H/D  =  0.28 
for  data  taken  at  Ri  >  5  and  H/D  =  0.3  for  Ri  >  1. 

(v)  Associated  with  the  outward- flowing  current,  there  is  an  inward  subsurface  jet  of 
velocity  U,  r  0.25  for  the  entire  Ri  range  investigated.  The  entrainment  flow  for  the 
expanding  plume  is  also  provided  by  this  jet. 

(vi)  When  Ri  <  20,  entrainment  of  heavy  fluid  from  the  lower  layer  into  the  plume 
impinging  on  the  interface  was  evident.  The  result  is  the  generation  of  a  sub-interfacial  horizontal 
current  directed  toward  the  region  of  impingement  The  velocity  of  this  current  was  found  to  be 
given  by  Ub/q^^^  =  0.2  Ri*l/^;  when  this  is  expressed  in  terms  of  a  vertical  entrainment  velocity 
Ug  pertinent  to  the  horizontal  area  of  impingement,  the  entrainment  law  Ug/Ool/^  -  o.6  Ri-*  could 
be  obtained.  Additional  entrainment  can  occur  into  the  gravity  current  when  Ri  <  1. 

(vii)  The  impingement  of  the  plume  on  the  density  interface  generates  interfacial  waves  of 
different  frequencies.  The  dominant  mode  was  found  to  be  the  lowest  iiKxle  of  internal  waves. 

The  above  experiments  were  extended  to  include  the  effects  of  background  rotation  by 
carrying  them  out  on  a  rotating  table.  It  was  demonstrated  that  the  leads-induced  motions  are 
affected  by  the  background  rotation  after  descending  to  a  depth  of  where  q^  is  the 

surface  buoyancy  flux  and  Q  is  the  rate  of  background  rotation.  The  width  of  the  plume  at  this 
point  is  1. 1  (qo/Q3)l/3.  After  some  time,  the  plume  becomes  unstable  at  its  transverse  edges  and 
deflects  sideways,  thereby  producing  a  strongly  three-dimensional  cyclonic  spiraling  flow  pattern. 
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